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Here we present the theoretical foundation of the strong coupling phenomenon between quantum 
emitters and propagating surface plasmons observed in two-dimensional metal surfaces. For that 
purpose, we develop a quantum framework that accounts for the coherent coupling between emitters 
and surface plasmons and incorporates the presence of dissipation and dephasing. Our formalism 
is able to reveal the key physical mechanisms that explain the reported phenomenology and also to 
determine the physical parameters that optimize the strong coupling. A discussion regarding the 
classical or quantum nature of this phenomenon is also presented. 



PACS numbers: 42.50.Nn, 73.20.Mf, 71.36.+C 

Surface plasmon polaritons (SPPs), hybrid bound 
modes comprising both electromagnetic fields and charge 
currents, are well-known to have both a subwavelength 
confinement and propagation lengths of tens or even hun- 
dreds of wavelengths [U [2]. For this reason, the inter- 
action between quantum emitters (QEs) and SPPs has 
attracted great interest recently 3 5 . It has been shown 
that QE-SPP coupling can lead to single SPP generation 
[BHS] and that the interaction between two QEs can be 
mediated by SPPs, resulting in energy transfer, super- 
radiance [9] and entanglement phenomena [T0HI2] . Re- 
cently, there have also been several experimental studies 
that show the emergence of strong coupling (SC), i.e., co- 
herent energy exchange between propagating SPPs and 
excitons either in organic molecules |13H18j or in quan- 
tum dots [19 21 . However, up to our knowledge, a first- 
principles explanation of these experimental results has 
not been presented yet. 

In this Letter we analyze the phenomenon of SC be- 
tween quantum emitters (or absorbers) and SPPs and 
present its theoretical foundation. We develop a com- 
plete quantum treatment that is able not only to calcu- 
late absorption spectra and reproduce the experimental 
phenomenology, but also to deal with more complex as- 
pects as photon statistics. 

In Figure 1 (a) we render a sketch of the general struc- 
ture that mimics the experimental configuration: a col- 
lection of N QEs immersed into a layer of thickness W, 
and placed on top of a thin metal film (thickness h). In 
this work the acronym QE will refer to a quantum sys- 
tem with discrete electronic levels, like organic molecules 
or quantum dots. In some of the experimental set-ups 
and in order to avoid quenching of the QEs, a dielectric 
spacer of width s is located between the QEs and the 
metal substrate. We will take e2 = ei = 1 in our cal- 
culations and we will use the dielectric function of the 
metal (silver), e m , as tabulated in Ref. 22J. As a mi- 
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FIG. 1: (color online) (a) Schematic picture of the N QEs 
distributed in a volume of width W separated by a distance 
s from a metal film of thickness h. (b) Coupling constant, 
gp(k; zo), for a single QE with perpendicular (solid blue) and 
parallel (red dashed) orientations (see inset), placed at zo = 
20nm and interacting with a SPP of momentum k(uj). 



nor simplification, we will assume a semi-infinite metal 
substrate instead of the metal film considered in the ex- 
periments (these films are thick enough for the SPPs to 
be very similar to those of a single interface). Each QE 
is represented by a two-level system (2LS) {\g) , |e)} and 
characterized by a transition frequency ujq (in this pa- 
per we will use ujq = 2cV, % = 1) and dipole moment /2 
with spontaneous decay rate 70 = WQ/i 2 /(37reoc 3 ). This 
description assumes a large separation between the elec- 
tronic energy levels of the emitter, with only one possible 
transition at the excitation frequency. 
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A QE placed in the vicinity of a metal surface can 
decay into three different channels [23]: excitation of 
SPPs that propagate along the metal surface, radiation of 
photons into the far field and dissipation through ohmic 
losses in the metal. In order to study the SC regime be- 
tween the QE and SPPs, the excitation of plasmons will 
be considered as the coherent channel, while radiation to 
free space and losses into the metal will be treated as 
dissipation mechanisms. The general Hamiltonian that 
describes the coherent interaction between N QEs and 
the 2D-SPPs can be written as (a detailed account of its 
derivation is presented in the Supplementary Material): 
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where aij and o~ \ ^ are the QE lowering and raising oper- 
ators of a QE that is located at (Ri, zj). In our modeling 
we assume the ensemble of QEs to be disposed in Nl 
layers, each of them having N s equal emitters such that 
N = Nl x N s . In Eq.(l), and at are the destruc- 
tion and creation operators for the SPP quantum field 
with in-plane momentum k and energy u>, linked by the 
dispersion relation k 2 (uj)(e m + 1) = e m ui /c . The area 
of the metal-dielectric interface is A and gp(k;zj) is the 
coupling constant of the dipolar interaction between a 
given QE and the SPP field: 
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where L(k) is the effective length of the mode [24j 125] . 
For calculating this coupling constant, propagation losses 
of the SPP modes are neglected. The unitary vectors in 
the k and z directions are ut and u z , respectively. The 
dependence of gp with z is dictated by the decay length of 
the SPP in the z-direction via k z = ^Jk 2 — ui 2 /c 2 . In Fig. 
1(b), we render the evolution of gp(k; z) with frequency 
for two possible orientations of the dipole: parallel to the 
momentum k and perpendicular to the metal surface. 
In both cases, the couplings are evaluated for QEs with 
7o = O.lmeV, which is a typical value for the J-aggregates 
used in the experiments as QEs [T3HT51 [17] . As shown in 
Fig. 1(b), the coupling constant between the QE and the 
SPP mode is larger for the perpendicular orientation, as 
k z is always smaller than |fc|. 

To simplify the general Hamiltonian (1), we first con- 
sider that, in the low excitation regime, the QE lowering 
and raising operators {pi,j and a\ ■), can be replaced by 

bosonic operators, bij and b\p respectively. Second, as 
in the experiments the ensemble of QEs is disordered, 



we assume that the structure factor is peaked at zero- 
momentum. Third, we build up a collective mode of the 
N QEs, Dt, by means of a transformation in which each 
excitation is weighted by its coupling to SPPs. Based 
on this, the total hamiltonian of the N QEs interacting 
with the SPP modes of a 2D-metal film can be written 
as H = J^k H% ( see t ne details of its derivation in the 
Supplementary Material), in which the hamiltonian as- 
sociated with momentum k has the following expression: 

£Tf = wol^s + ^(fc)4«£ + \9$ $)H D \ + c - c -l- ( 3 ) 
Here, g^ (k) is the effective coupling constant: 
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The last equality in Eq. (4) assumes a continuum of 
layers in the z-direction with a total thickness W and a 
volume density of emitters n = N S N l / [AW) . The hamil- 
tonian as written in Eq. (3) is one of the main results 
of our work, as it allows an ab-initio quantum treatment 
of the coherent coupling between an ensemble of N QEs 
and SPPs. Notice that this interaction conserves the to- 
tal momentum of the system composed of the supermode 
of QEs and the SPP. When evaluating the coupling con- 
stant for a momentum k, gj~ (k), there is no need to rely 
on fitting parameters and can be calculated from first 
principles, as shown below. 

In Figure 2(a) we plot the effective coupling con- 
stant g N evaluated at k(u>o) for a density of emitters 
n = 10 6 /im~ 3 (of the order of the volume densities used 
in the experiments) as a function of W and for differ- 
ent values of the spacer width. This magnitude de- 
pends on the orientation of the QEs' dipole moments. 
Here we render the two limiting cases (all dipoles ori- 
ented perpendicularly or parallel to the metal surface) as 
well as an isotropic average over these two orientations, 
g 2 so = 2g 2 /3 + g\/i. Two main conclusions can be ex- 
tracted from this figure. First, g N depends strongly on 
W but saturates for thick enough films. This saturation 
is due to the exponential dependence of gp on z, related 
to the spatial decay of the SPP mode, and, therefore, 
is determined by the dielectric environment of the metal 
film. Second, the dependence of g^ (k) on the width of 
the spacer layer is not very strong. 

The excitation of the hybrid system needs to be in- 
cluded into the theoretical framework. In order to re- 
produce the typical experimental configuration, we will 
assume that SPPs are excited by a coherent laser field. 
A new term is incorporated into the total Hamiltonian, 
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sures the intensity of the laser field and is the oper- 
ating frequency of the laser. In this way, the laser field 
fixes the SPP parallel momentum, k, implying that only 
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ciated with each of the three dissipative channels. Re- 
calling that the general expression of a Lindblad term 
associated with an arbitrary operator c is C c — (2cpc' — 



c^cp 



pc^c), the master equation for the density matrix 



associated with momentum k, pt{t), can be written as: 
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The solution of the master equation for fco = k(uo) 
(in-plane momentum that displays maximum coupling) 
yields to coherence functions being proportional to 
exp(iitt) where R is the Rabi splitting at resonance: 



R = 
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FIG. 2: (color online), (a) Coupling constant, g (W), for sep- 
arations s ranging from 1 to 50 nm and for parallel, perpendic- 
ular and isotropic orientations of the QEs with 70 =0.1 meV. 
(b) Real (solid black) and imaginary (red dashed) parts of 
the Rabi splitting at resonance for dipoles oriented isotropi- 
cally, Ri S o, as a function of n for the geometrical parameters: 
s — lnm and W = 500 nm and 7^ = 40meV. (c) Polariton 
population (see main text) of a distribution of QEs as a func- 
tion of k, with the same geometrical parameters as in panel 
(b) and with Og — O.lg . The volume density in this case is 
n = 10 6 pm~ 3 , as in panel (a). 
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taken into account. 

Finally, the description of the dynamics of the sys- 
tem must be completed by considering both the losses in 
the ensemble of QEs and the dissipation associated with 
the SPP mode. The decay lifetime of the SPP mode, 
7 a _, can be calculated from the SPP propagation length, 



Following the standard analysis [55], we will consider 
that our hybrid system is within the SC regime when the 
imaginary part of the Rabi splitting is zero. In Fig. 2(b), 
we plot the evolution of R = R r + 1R4 with the volume 
density n for an ensemble of N QEs whose dipoles are 
oriented isotropically. For very low densities (for this set 
of parameters, n < 2 x 10 3 /im -3 ), R is a purely imagi- 
nary number and, therefore, the system operates in the 
WC regime. This density threshold, n tl is mainly con- 
trolled by 70 as 7^ 3> id, la for this set of decay rates. 
Notice that as 7^ decreases exponentially when lowering 
the temperature [2B], n t is expected to be much smaller 
at very low temperatures (by assuming 7^ = at zero 
temperature, n t would be around 20 pm~ 3 ). For high 
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typical densities 



in the experiments [18]), R r (the so-called vacuum Rabi 
splitting) is dominated by the coupling constant g N as 



Lgpp, and group velocity, 
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lifetime increases as the frequency approaches the SPP 
cut-off frequency, being around 5 meV for uj = luq = 2 eV. 
The lifetime associated with the collective mode of the 
ensemble of N QEs, 7d £ , is obtained from the averaged 
value of the decay rates for each individual QE, 7 CT (z), 
weighted by a term proportional to \gp(h, z)\ 2 (for details 
see Supplementary Material). Additionally, in order to 
be as close as possible to the experimental conditions, the 
existence of vibro-rotational states in organic molecules 
must also be taken into account. These degrees of free- 
dom within the QEs can be incorporated into the 2LS 
model by means of pure dephasing mechanisms, charac- 
terized by a dephasing rate, 7^. In this work we take 
70 = 40 meV |26j , which is a typical value at room tem- 
perature for the organic molecules used to observe SC 
between N QEs and SPPs. 

With all these ingredients, we use a Markovian master 
equation for the density matrix and introduce pertur- 
batively the corresponding Lindblad operators [57] asso- 



g" {7_D,7aj70j an d Rr ~ 9 ■ As this coupling con- 
stant scales as y/n, so does R r , as observed in the ex- 
periments. Within our formalism, it is also possible to 
evaluate the absorption spectra, a magnitude that is at- 
tainable experimentally. In Fig. 2(c) we plot the polari- 
ton population (the sum of both the QEs supermode and 
SPP mode occupations, a magnitude that is proportional 
to the absorption by the system [5S]) versus energy and 
parallel momentum, showing the anti-crossing between 
the flat band at wo associated with the collective mode 
of the N QEs and the dispersive band of the SPPs. Al- 
ready existent experimental results [T5] can be confronted 
with our theoretical framework. In that experiment, the 
metal film was silver, W = 50nm and an ensemble of 
n = 1.2 x 10 8 /im~ 3 Rodhamine 6G molecules were used 
as QEs (70 = IpeV, as reported in Ref. [30]). This re- 
sulted in the observation of a Rabi splitting of 0.115 eV. 
For those parameters, our theory predicts R r — 0.04 eV 
for parallel-oriented QEs, R r — 0.18 eV for the perpen- 
dicular orientation and R r = 0.10 eV for an isotropic 
average, showing a good agreement between theory and 
experiment. 

Now we address the fundamental question regard- 
ing the classical/quantum nature of the SC regime ob- 
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served in this type of systems. Although a semi- 
classical formalism fed with phenomenological param- 
eters is able to reproduce qualitatively the reported 
absorption spectra [3T], this should not be taken as 
a statement that the system contains no interesting 
quantum physics. Non-classicality is unambiguosly 
revealed by the presence of photon antibunching in 
the dynamics of the strongly coupled system. For 
this reason, we analyze the behavior of the second- 
order correlation function, g^ 2 \ defined as g^{r) — 
lim t ^(Dl(t)(DlD s )(t + T)Dz(t))/(DlD~ k (t)y. Photon 

antibunching yields <7^(0) < 1. However, within the 
approximations leading to Hamiltonian (3) with con- 
structed from bosonic operators, the system behaves as 
two coupled harmonic oscillators. In this case, g^ (0) is 
always greater or equal to 1 [35] and its time evolution 
critically depends on the excitation means. For the case 
of coherent pumping, the system acquires the statistics 
of the laser field and hence g^(r) = 1. The case of inco- 
herent pumping is simulated in our theoretical formalism 
by introducing a Lindblad term, Pu k C D ^ [2H]) into the 

k 

master equation [Eq.(5)], instead of Hi. As shown in the 
inset of Fig. 3(b) (green-dashed line), when the collective 
mode is driven incoherently </ 2 ' (0) = 2 and its time de- 
pendence presents some Rabi oscillations but g^ 2 \r) > 1. 

In order to find fingerprints of non-classicality in our 
system (g^ 2 '(0) < 1), it is necessary to incorporate a non- 
linear term into Hamiltonian (3): 

H nl = JZUnDljA j?^ . (7) 

k.k' ,q 

The physical origin of this term can be twofold |33j : a 
direct coupling between the QEs, similar to the coulomb 
interaction between excitons reported in semiconductor 
structures, and/or saturation effects. In this last case we 
can even quantify this contribution by introducing the 
second-order correction within the Holstein-Primakoff 
approach [31] in the process of replacing the QE lowering 
and raising operators by the bosonic ones. By considering 
a quasi-2D layer of QEs, a non-linear term as expressed 
by Eq. (7) can be straightforwardly obtained from the 
general Hamiltonian (1) with Ud being —ojq/N (tech- 
nical details are given in the Supplementary Material). 
Notice that whereas in the linear case the key parameter 
is the volume density n, the saturation contribution to 
the non-linear term is controlled by the total number of 
active QEs, N. In Fig. 3(a), we show the dependence of 
<?^(0) on \Ud\ and the frequency detuning, {loq — u)l), 
both expressed in units of g N . In these calculations we 
have taken a dephasing rate 7^ = in order to find the 
most favorable, yet experimentally feasible, conditions to 
observe photon antibunching. As we consider pumping 
to only one fc-state, the population of a SPP-mode with 
parallel momentum k' is proportional to 8{k~k'), cancel- 
ing out the summation in k' in Eq.(7). In addition, the 



summation in q can also be neglected because the shape 
of the SPP dispersion relation does not allow parametric 
scattering j35j [36] , in which both energy and momentum 
are conserved, to SPP states with q ^ 0. Two partic- 
ular cases (|!7d| = 0.025^ and 0.005^^) are displayed 
in Fig. 3(b) for a better visualization. If we assumed 
a saturation origin for Ud, these two cases would corre- 
spond to N R3 2 x 10 3 and N ~ 10 4 , respectively. Im- 
portantly, photon antibunching is observed in both cases 
and is greater when the laser frequency almost coincides 
with u>o or uiQ±g N . Therefore, our results suggest that in 
order to observe noticeable photon antibunching, the ex- 
periments should be performed at very low temperature 
to avoid dephasing and plasmon losses. Additionally, in 
order to reduce the number of active QEs, the laser beam 
should have a very small spot size, and the QEs should 
be disposed forming quasi-2D layers. 
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FIG. 3: (color online) (a) Contour plot of JT 2 '(0) as a function 
of the non-linearity, \Uo\/g N , and detuning, (wo — uJh)/g N , 
for the coherently pumped configuration, with system param- 
eters: s = 10 nm, W = 10 nm and n = 10 6 /im -3 , which 
yields g N w 50 rneV, see Fig. 2(a). The color code is: blue; 
1 white; 2 red. (b) Horizontal cuts of (a) at two fixed non- 
linear parameters: \Ud \ — 0.005(7^ (dashed red) and 0.025(7^ 
(solid black). Inset: g {2) (r) for a system with \U D \ = Q.Q2hg N 
and wl = wo +Q.lg (solid black) and with |J7d| = for both 
incoherent pumping (dashed green) and coherent excitation 
(dotted blue). In all these calculations we have taken a de- 
phasing rate 7^ = 0. 

In conclusion, we have presented an ab-initio quantum 
formalism to study the phenomenon of SC between QEs 
and propagating SPPs in two-dimensional metal surfaces. 
Based on this formalism, we are able to predict the crit- 
ical density where SC emerges for a given geometry and 
distribution of QEs, and to determine the optimal geo- 
metrical parameters that maximize SC. Our results show 
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that, for experiments carried out at room temperature, 
QE and SPP losses play a minor role in the emergence of 
SC. Both coherent coupling between the QEs and SPPs 
and pure dephasing mechanisms determine the strength 
of the phenomenon in this case. Additionally, the devel- 
opment of this general quantum framework allows study- 
ing the fundamental nature (classical versus quantum) of 
this phenomenon by analyzing the conditions in which 
photon antibunching could be observed. 
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